Solutions of the Gaudin Equation and Gaudin 

Algebras 

A. B. Balantekin* 
Department of Physics, University of Wisconsin 
Madison, Wisconsin 53706 USA 

T. Derelit 

Department of Physics, Koq University 
34450 Sariyer, Istanbul, Turkey 

Y. Pehlivan * 

Department of Mathematics , Izmir University of Economics 
35330 Izmir, Turkey 

11 March 2005 



Abstract 

Three well-known solutions of the Gaudin equation are obtained under 
a set of standard assumptions. By relaxing one of these assumptions 
we introduce a class of mutually commuting Hamiltonians based on a 
different solution of the Gaudin equation. Application of the algebraic 
Bethe ansatz technique to diagonalize these Hamiltonians reveals a 
new infinite dimensional complex Lie algebra. 
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1 Introduction 



The Bardeen-Cooper-Schrieffer (BCS) pairing model Hamiltonian was diag- 
onalized by Richardson in 1963 in an algebraic way The step operators 
introduced in Richardson's solution belong to an infinite dimensional com- 
plex Lie algebra which is usually referred to as the rational Gaudin algebra. 
It is one of the three algebras which emerged from Gaudin 's work during 
the 1970's [2] who introduced a new class of integrable models starting with 
the following operators: 

TV 2 



where wfj are complex numbers to be determined and tf are the generators 
of N, mutually commuting SU (2) algebras. The latter obey the standard 
SU(2) commutation relations: 

[t+,tj] = 2S ij t° j , [tltf] = ±6 i:j tf , i,j = 1,2,. ..,N (2) 

where tf = t\ ± itf. Gaudin showed that the operators in Q mutually 
commute 

[hi,hj] = (3) 
if and only if the coefficients satisfy the equations: 

<™]k + < w l ~ = ( 4 ) 

for all distinct triples (i,j,k) and for all permutations of the upper indices 
(0,1,2). There are three solutions to equation (j3J) under the following as- 
sumptions: i) The coefficients wfj are antisymmetric under the exchange of 
the indices i and j 

«,?. + u# = 0, (5) 

ii) each coefficient wfj can be expressed as a function of the difference be- 
tween two real parameters and Uj, and iii) each hi commutes with the 
operator T = Y^j tj so that the z component of the total SU (2) is conserved. 
Gaudin's solutions are given as follows: 

v& = for a = 0, 1, 2, (6) 

J Uj — Uj 



w%=pcot[p{u l -u J )} , wl j =wf j = -r—^ TTvf ' ( 7 ) 
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«4=P«rth|p(ii i -u f )] , < = < = sinm l_ Uj) y (8) 

which are commonly referred to as rational, trigonometric and hyperbolic 
solutions, respectively. Here, p is a real parameter . 

The operator which is obtained by substituting the rational solution into 
(PQ) is called rational Gaudin magnet Hamiltonian: 



N 



(r) _ t i ■ t j 



3=1 



where 



t i ■ t j — tjtj + t i tj + t j tj — tj + + i j 



2 1 3 1 3 ' 

Similarly, substituting the trigonometric and hyperbolic solutions into 
we obtain 



N 



+ ( 10 ) 



5=1 



N 



2sin[p(uj - Uj)] 



3=1 



which are called trigonometric and hyperbolic Gaudin magnet Hamiltonians, 
respectively. 

Application of the algebraic Bethe ansatz method to the Gaudin mag- 
nets yields the rational, trigonometric and hyperbolic Gaudin algebras E] • 
These are infinite dimensional complex Lie algebras which are related to par- 
ticular solutions of the classical Yang-Baxter equation 5 . As a result, each 
algebra admits a one-parameter family of mutually commuting Hamiltoni- 
ans H(X). Here A is a complex parameter which is usually referred to as the 
spectral parameter. These operators may be identified as the integrals of 
motion of a quantum system as well as the traces of transfer matrices of a 

1 These solutions are not completely unrelated. Gaudin equation Q is satisfied for all 
complex values of p. But only for the real and pure imaginary values of p Hamiltonians 
in Q are Hermitian. On the other hand, substituting ip in place of p, one can convert 
trigonometric and hyperbolic solutions into each other. For this reason, we restrict our- 
selves with the real values of p. Also note that in the limit p — ► trigonometric and 
hyperbolic solutions go to the rational solution. 
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vertex model. Main objective in both cases is to diagonalize them simulta- 
neously and this is achieved by either functional or algebraic Bethe ansatz 
techniques 6 . Richardson-Gaudin methods have recently found many ap- 
plications in quantum many-body physics |7j. 



2 Further Solutions of the Gaudin Equation 

In this paper we investigate a different solution of the Gaudin equation Q . 
Under the three assumptions listed in the introduction, all possible solutions 
are enumerated above. To find new solutions we need to relax one or more 
of these assumptions. We keep the constraints (ii) and (hi) but generalize 
the constraint (i) to 

< + w % = (12) 
where q is a real parameter 2 . This solution is given by 

wfj = q coth[q(ui — Uj)] — q for a = 0,1,2. (13) 

The operators which are obtained by substituting this new solution in Q 
will be denoted by lif ' ': 

N 

= Yl fa - Uj )} -q)~ti- ~~tj. (14) 

3=1 

Since ((T3"|) is a solution of Gaudin equation (JI|, these operators mutually 
commute: 

[h?,h < f ) ] = 0, (15) 

In the limit q — > 0, the new solution given by (|13|) approaches Gaudin's 
rational solution Q. As a result, 

lim/if =h i [ ) . (16) 
g->0 

The solution presented here is a different solution of Gaudin equation 
and the Hamiltonians given by (|14|) for q ^ cannot be obtained from the 
rational, trigonometric or hyperbolic Gaudin magnet Hamiltonians. To see 
that, let us assume for a moment that there exists an operator S such that 

Sh^S- 1 = h[ q) (17) 



2 In general, Jl} is satisfied for all complex values of q. But only for real values of q one 
obtains Hermitian Hamiltonians in ijljl. 
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where h\ represents any one of the Hamiltonians given in equations (jyj)- 
Suppose we take the sum of both sides of (|17j) : 



N N 

,(*)<y-l _ h(i) 



J2Sh^S-i=Y,hf. (18) 



i=l i=l 

Antisymmetry of © implies 

N 



E^ ) = ( 19 ) 

i=l 



for rational, trigonometric and hyperbolic magnet Hamiltonians. Hence the 
sum on the left hand side is equal to zero. On the other hand, the sum on 
the right hand side of (|18|) is not equal to zero. Instead the new condition 
given in (|12ft implies 



N N 



E h i = -9 E *<• *i« ( 20 ) 

i.e. an operator 5 satisfying ((T7|) does not exist for q ^ 0. 



3 Gaudin Algebras 

If one applies the algebraic Bethe ansatz technique to diagonalize the Gaudin 
magnet Hamiltonians which correspond to the solutions © , , or one 
uses the step operators belonging to the rational, trigonometric or hyper- 
bolic Gaudin algebras, respectively. Below we will show that if the algebraic 
Bethe ansatz method is employed to diagonalize the Hamiltonian of (|14|). 
then a different algebra emerges in a completely analogous way. Before in- 
troducing this algebra, however, we would like to outline the basic features 
of rational Gaudin algebra and its relation to the rational Gaudin magnet 
Hamiltonian, as an illustration of the technique. Trigonometric and hyper- 
bolic Gaudin algebras and their relation to the trigonometric and hyperbolic 
Gaudin magnet Hamiltonians follow very similar lines. 

Rational Gaudin algebra is generated by three families of operators 
J + (A), J~(A) and J°(A) parameterized by a complex number A. Commuta- 
tors are given by 

A — fj, 
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[J°(X),J^)] = ± J±iX) ~ J±(t " ) , (21) 

A — fx 

[J°(X),J°^)] = [J ± (X),J ± ^)]=0. 

Commutators of J + (A), J - (A) and J°(A) at the same value of the complex 
parameter are given by taking the limit li — > A. From these commutation 
relations it is easy to show that 

H(X) = J°(A)J°(A) + Ij+(A)J-(A) + ij~(A)J+(A) (22) 

form a one-parameter family of mutually commuting operators: 

[H(X),H(li)]=0. (23) 

Starting from a lowest weight vector and using J + (A) as step operators, one 
can diagonalize these operators simultaneously. Lowest weight vector |0 > 
by definition satisfies 

J-(A)|0>=0, and J°(A)|0 >= W{X)\0 > (24) 

for every A. Here W(A) is a complex-valued function. The state |0 > itself 
is an eigenvector of H(X): 

H(X)\0>= £ (A)|0> (25) 

with the eigenvalue 

E (X) = W{Xf - W'{X) (26) 
where the prime denotes derivative with respect to A. In addition 

16,6, ■■■^n>= J + (6)^ + (6) • • • </ + (60|o > (27) 

is an eigenvector of H(X) if the quantities 6, 6, •••,6 satisfy the set of 
Bethe ansatz equations 

n 

^«)=Ei r for « = l»2,...,n. (28) 

J^i &» - 4/3 

(For an analytic solution of these equations see e.g. References 8 and 9 ). 
The corresponding eigenvalues will be 

W^W-^ MM . (29) 
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There exists a realization of the rational Gaudin algebra in terms of the 
SU(2) generators of ©, given by 



r-f «i - A f-f «j - A 

Here u\,U2, ■ ■ ■ ,un are arbitrary real numbers which are all different from 
each other and N is a nonnegative integer 3 . Assuming that the eigenvalues 
of the Casimir operator of the j th SU(2) are Sj(sj + 1), the corresponding 
W(X) is given by 

N 

W(A) = £^\. (31) 
In this realization H{\) is given by 

H(\) = V -. -. (32) 

We see that H(\) has simple poles on the real axis. Residues of H(\) at the 
points A = Ui are proportional to the rational Gaudin magnet Hamiltonians 
given in ©: 

-^Res{i?(A))} A=Ui =/4 r) (33) 

Equation (|2*3^1 implies 

[H(X),h ( f ) } = 0, and [hf\hf]=Q. (34) 

As a result, rational Gaudin magnet Hamiltonians h\ are automatically 
diagonalized together with H{\). It is easy to read off the eigenvalues of 
the rational Gaudin magnet Hamiltonians from the eigenvalues of H(X) as 
follows: 

^ (r ' ) |0>=4o ) l°>' (35) 

where 

N 

4? = £ ( 36 ) 

^ - Uj 

J=l J 

3 In general iti,lt2,.-. ,wjv can be complex and they need not be different from each 
other. But in most physical applications we are interested in those realizations for which 
ui,U2, . . . ,un are real and different from each other. Reality guarantees that J + (A)^ = 
J" (A*) and J°(A)t = J°(A*) 
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and 

h { [ ] 16, 6, ■•.,£«>= ■ ■ ■ > &» > ( 3? ) 

where 

a=l Ul 

Here |0 > is the lowest weight state and £ a are the solutions of the Bethe 
ansatz equation (|28[). 

We next search for an operator H^ q \X), the residue of which gives us 
the operators nf (cf. Equation (J33J)). We wish to write such an operator 
using an algebra similar to that in (|21j) . Below we show that such an algebra 
exists and its generators J q (X), J q (X) and J q (X) satisfy the commutation 
relations given as 

W"(A), JJ-O*)] = 2 g £hL(A J - ( ^)] + 29 (J ° (A) + J > ]) ' 

[JJ(A), Jf (M)] = ± ^fl A ^ A J !^ ] ) ± ? WW + (M)) , (39) 

[J 9 (A),J0( M )] = [J±(A),J 9 ± (^)]=0. 

We observe that, in the limit q — > 0, commutators in (|39|) approach the 
commutators of Gaudin algebra given by (|21j) . It should nevertheless be 
emphasized that the algebra of (|39|) is not a q-deformed algebra, but an 
ordinary Lie algebra. We can show that 

tfW(A) = J° q (X)J° q (X) + ^J+(A)J-(A) + ^J-(A)J+(A) (40) 

form a one-parameter commutative family: 

[H^(X),H^(fi)] = 0. (41) 

We can simultaneously diagonalize them starting from a lowest weight vector 
satisfying 

J-(A)|0>=0, J?(A)|0>=W ff (A)|0> (42) 

where W q (X) is a complex valued function. It can be shown that |0 > is an 
eigenvector of H^(X) with the eigenvalue 

4 9) (A) = W q (X) 2 - WAX) - 2qW q (X). (43) 
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In general, one can write a Bethe ansatz form for the eigenvectors of H^ q > (A) 
as 

16,6, ...,Zn>= (6) ■ ■ ■ 4(^)\0 > ■ (44) 

For this to be an eigenvector of H^ g >(X), the complex numbers 6, 6, • • • , 6 
must satisfy the following system of equations: 



W q (Ca)= E ^M* 1 [?(£«" &)]-!) for a = l,2,...,n. (45) 

/3=1 

In this case, eigenvalue associated with the state 6, • • • > 6 > is given by 



4«)(A) = 4 9) (A) - 2 £ g (coth [q(X - 6)] - 1) (W q {\) - W q (U) . (46) 

a=l 



This algebra admits two realizations in terms of the SU(2) generators of (J2J). 
These are given by 



N 

1=1 

and 

iV 

J P°w = E « ( tanh - A )i + !) ^ ,0 - ( 48 ) 

i=l 

In the limit g — > 0, the first realization goes to the realization of the rational 
Gaudin algebra given by ()3Uj) whereas the second realization vanishes. 
In the first realization, Equation ((171) . the operator H^ q \X) becomes 

N 

H {q) (A) = E ^ coth ^ ~ A )] + 9) (« coth fe( u i - A )] + l) *i ' *j • (49) 

We see that H^(X) has simple poles on real axis at the points A = itj. It is 
easy to show that —1/2 times the residue of H^ q \X) at the point A = Ui is 

-^Res {# ? (A))} A=Ui = ^ - qt l ■ ti. (50) 

Since tj ■ ti commutes with every hf and H^(X), (JUj) implies 

[#(«)(A),/4 9) ] =0, and = 0. (51) 



S 



As a result, the operators fif^ are automatically diagonalized together with 
i?(«)(A). As in the case of the rational Gaudin algebra, one can read off the 
eigenvalues of hft from the eigenvalues of H^(X): 



I$>\0>=E$\0> t (52) 



where 

and 
where 



N 



E i = E s i s jl(. coth - [<l( U i ~ u j)} ~ !) ( 53 ) 



hf ] |6, 6, ■ •-,£«>= 16, 6, • • • , 6 > (54) 

n 

= E$ ~ s * E «( coth - &)] - !)■ ( 55 ) 

a=l 

Details of the derivation of these eigenvalues are given in the Appendix. 

4 Linear r-Matrix Structure 

The rational, trigonometric and hyperbolic Gaudin algebras have a linear 
r-matrix structure. This means that one can write down the commutators 
of these algebras in the following matrix form: 

[L(A) <g> I, I <g) L(ji)] + [r(A - /i), L(A) (g)I + I<g) L(/x)] = (56) 

where / is the 2x2 identity matrix, r(A — fi) is the r-matrix described below 
and L(X) is given by 

J°(A), J + (A) and J - (A) are generators of the rational, trigonometric or hy- 
perbolic Gaudin algebras. In ifHHj) . the term [Z(A)®I, is equal to a 
4x4 matrix whose elements are various commutators between the generators 
J°(A), J+(A) and J - (A). On the other hand r(A-/i) is a Mat 2 (C)(g)Mat 2 (C) 
matrix valued function of A—// which carries the information about the struc- 
ture constants of the algebra. For instance, for the rational Gaudin algebra, 
the r-matrix is given by 

r(A -/*) = —!— P (58) 

A — fl 
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where P is the permutation matrix on 



which is equal to 



P 



It is customary to introduce the following three mappings 



A 
















1 








1 








Vo 











(59) 



(60) 
(61) 
(62) 



where a, b are 2x2 matrices and I is the 2x2 identity matrix. Then r y 
are defined as follows: 



V 12 


a®b — 


-* a ^ 


56$ 


5/, 


if 13 


a®b — 


-> a £ 






cp 23 


: a<g>b — 




<D a C 


3 b, 



,12 



.13 



¥> id [r(A-«r)], 



„23 



^[r(/i-<r)]. (63) 



Since r-matrix carries information about the structure constants of the al- 
gebra, Jacobi identity leads to the equality 



[r 13 ,r 23 ] + [r 12 ,r i3 ] + [r 



12 „13n 



42 231 







(64) 



which is the well known classical Yang-Baxter equation (for solutions of this 
equation see [TU|, for a review see ^U). Equation (|64|) guarantees the 
mutual commutativity of the operators H(X) which can now be written as 



H(\) = ±Tr{L(\f} 



(65) 



Trigonometric and hyperbolic Gaudin algebras also have the r-matrix struc- 
ture described above and their r-matrices are also solutions of the classical 
Yang-Baxter equation. 

In order to see if one can repeat the same procedure for the algebra we 
introduced in (39), we study a matrix of the form 



L q (X) 



fJ° q (X) J+(A) 



(66) 



\J+(X) -JJ(A), 
We then substitute it into the equation 

[L q (X) ® I, I ® Lgifj,)] + [r(A - /i), L q (\) ® I + 1 ® L,(p)] = 0, (67) 
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assuming the existence of such an r-matrix. One can compute the left-hand 
side of this equation and then set all the components equal to zero to find the 
components r tJ of this r-matrix. For instance, the (21) and (13) elements 
are equal to 

0"22 - rn + gcoth [q (A - li)] - q) J q (li) 
+ (r 23 - q coth [q (A - li)] - q) J- (A) - r 4 i J+ (A) + 2r 2i J° (li) 

and 

(rn ~r 3 3-q coth [q (A - Li)] - q) J+ (A) 
- (r23 - gcoth [q (A - li)] + q) J+ (li) + r w J 9 " (/i) - 2r i3 J° (A) , 

respectively. Setting these elements equal to zero requires setting the coeffi- 
cients of all the algebra elements equal to zero. The coefficient of J q (\) in 
the first equation and the coefficient of J q (li) in the second equation gives 

r 2 3 - q coth[q(A — //)] — q = and r 23 - q coth[g(A — li)] + q = 0, (68) 

respectively. But these equations are incompatible unless q = 0. Repeating 
the same calculations for the other components, one can conclude that there 
is no r-matrix which satisfies (|67|) . Instead of the L q (\) matrix in (|66j) . one 
may try to find a more general form of L q (\) matrix for which an r-matrix 
can be found to satisfy (|67|). But we were unable to find such an L q (X) 
matrix 4 . On the other hand we were able to show that one can write the 
commutators of the algebra introduced in ()39|) in the following form: 

[L g (A)®J,I®i g (M)]+[r g (A-M),ig(A)®i] + [r_g(A-/i),/®L g (At)] = 0. (69) 

Here L q (\) is given by (|66|) and the r^-matrix is given by 

r,(A -Li) = (q coth [q(X - li)] + q)P. (70) 

P is the permutation matrix introduced in (|59jl. As in the case of the 
r-matrices of rational, trigonometric and hyperbolic Gaudin algebras, the 
r g -matrix carries information about the structure constants of the algebra 

4 For instance, for a matrix of the form 
L q {\) 



a, j 



J°(A) /3J+(A)+ 7 J-A) 
PJ~{\) -7J + (A) -aJ°(X) 

Tr[L q (X)Lq(\)] is proportional to H^- q \X) when a 2 = /3 2 + 7 2 . However, one can similarly 
show that there is no solution to I67|l for this L q {\) either. 
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of 1)3*9")) . r % q matrices, defined by similar maps as in Equations l|60 |) -l|63 |) . 
satisfy the following equation: 

V% ff] + \r% rf] + [r% rf] = 0. (71) 

We see that the r-matrix structure of the algebra introduced in (|39|) is not the 
same as the r-matrix structures of the rational, trigonometric and hyperbolic 
Gaudin algebras. Nevertheless, it shares with the other three algebras, the 
crucial property of admitting a one-parameter family of mutually commuting 
operators iJ^(A). Similar to 1)65(1 . one can write this operator as the trace 
of the square of L q (X): 

H^{\)= l -Tr{L q {\?}. (72) 

Note that in the limit q — > 0, the algebra of 1)39)1 approaches the rational 
Gaudin algebra. In this limit both r q and r_ 9 go to the r-matrix of the 
rational Gaudin algebra given in ()58|) and from 1)71)1 we recover the classical 
Yang-Baxter equation (|fi4jh 



5 Conclusion 

In this article, by relaxing one of the conditions imposed by Gaudin, we 
presented a different solution to the Gaudin equation and wrote down the 
corresponding set of mutually commuting Hamiltonians. We also identified 
the related infinite dimensional Lie algebra. This algebra allows a one pa- 
rameter family of mutually commuting Hamiltonians H^ q ' (A) parameterized 
by a complex spectral parameter. We diagonalized these Hamiltonians us- 
ing the algebraic Bethe ansatz method and constructed the eigenvectors and 
eigenvalues. The procedure is parallel to the other known Gaudin algebras. 

A different generalization of the Gaudin algebras was recently introduced 
in Ref. |13j that is based on the rational, trigonometric, and hyperbolic so- 
lutions of the Gaudin equations. Our generalization is based on a different 
solution of the Gaudin equation obtained under different assumptions. 
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Appendix: Computation of the Eigenvalues 

To see that the lowest weight state given by (|42jl is an eig enstate ofH^(X) 
we first use the commutators (|39|) to write H^(X) in the following form: 

Ft«>(A) = J°(A)J°(A) + J + (A)J-(A) + i[J'(A),J + (A)] 

= J°(A)J°(A) + J+(A)J-(A) (73) 

_I i im 2 q ( J ° ( ^ ~/° {fI l + J°(A) + J\n)) ■ 
2^\ H \tanh[q(X- fi)} y ' J 

Then the action of H^- q \X) on the lowest weight state is 

ffW(A)|0 >= (W(\) 2 - W'{\) - 2qW{\)) |0 > . (74) 

We see that the lowest weight state is an eig enstate of H (q \X) with the 
eigenvalue Eq{X) given by 

Now consider the state |£ >= J + (£)|0 >• Action of H W (A) on J+(£)|0 > 

is 

ff(?)(A)J+(O|0 >= 2q (coth [g(A - £)] - 1) W(£) J+(A)|0 > (75) 
+ (£ (A) - 2gW(A) coth [q(X - £)] + 2qW(X)) J+(£)|0 > ■ 

We see that H^(X) J + (^)|0 > is a superposition of the states J + (A)|0 > 
and J + (£)|0 >. Therefore, for a generic £, the vector J + (£)|0 > is not an 
eigenstate of H^(X). But if £ is a root of the TV (A) then the coefficient of 
J + (A)|0 > vanishes. In this case J + (£)|0 > is an eigenstate of H^ q \X) with 
the eigenvalue Ei(X) given by 

E[ q) (X) = E (X) - 2 (gcoth [q(X - £)] ~ q) W(X). (76) 
This is the energy one finds by substituting n = 1 and W{^) = in Eqn. 

For n > 1, action of H^(X) on the state J+(6)^ + (6) • • • J + (£n)\0 > is 
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given by 



ff<«)(A) J+(6) • • • «/ + (en)|o >= J + (6) J+(£ 2 ) • • • J + (£„)# (9) (A)|0 > 



2<7 (coth [q(X - &)] ~ 1) 



a=l 



\ 



Y,Q (coth [g(£ a - &)] - 1) - Wfa) 



J 



X J+fa) . . . J + (^-l)J+(A)J+(^ +1 ) . . . J + (£n)|0 > 



+ ^2g(coth [ g (A -&)]-!) 



a=l 



\ 



Xycoth[g(e«-&)]-i)-w(A) 

/3=1 



Xj+(ei)J + (6)... J + (Cn)!0> 



(77) 



Here we see that i?w(A) J + (£i) J + (^) • • • <^ + (£n)|0 > is a superposition of 
the states 

J + (6) • • • J+(£ a -i) J + (A) J+(£ a+ i) • • • J + (£n)|0 > (78) 
for a = 1, 2, . . . , n and the state 

j + (ei)j + (6)-..«/ + (en)|o> . (79) 

Coefficients of the states J + (£i) . . . J + (A) . . . J + (£ n )|0 > vanish when the 
conditions are satisfied. Consequently, J + (£i) J + (^2) • • • J + (Cn)\0 > 
become an eigenstate of H ( - q \\) and its eigenvalue is given by (|46|). 
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